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Abstract. The zero temperature phase diagram of binary boson-fermion mixtures 
in two-colour superlattices is investigated. The eigenvalue problem associated with the 
Bose-Fermi-Hubbard Hamiltonian is solved using an exact numerical diagonalization 
technique, supplemented by an adaptive basis truncation scheme. The physically 
motivated basis truncation allows to access larger systems in a fully controlled and very 
flexible framework. Several experimentally relevant observables, such as the matter- 
wave interference pattern and the condensate fraction, are investigated in order to 
explore the rich phase diagram. At symmetric half filling a phase similar to the Mott- 
insulating phase in a commensurate purely bosonic system is identified and an analogy 
to recent experiments is pointed out. Furthermore a phase of complete localization of 
the bosonic species generated by the repulsive boson-fermion interaction is identified. 
These localized condensates are of a different nature than the genuine Bose-Einstein 
condensates in optical lattices. 
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1. Introduction 

Ultracold atomic gases in optical lattices have become a powerful and flexible 
experimental tool for studying strongly correlated quantum systems. The unique degree 
of experimental control makes them an ideal model system to investigate fundamental 
quantum phenomena which have an impact on several fields, ranging from condensed 
matter physics to quantum information processing. Theoretically these systems are well 
described in the language of the Hubbard model. For bosonic gases the Bose-Hubbard 
model [T] led to the prediction of a superfluid to Mott-insulator transition in ultracold 
atomic gases by Jaksch et al. |2j. Since the experimental veriflcation of this quantum 
phase transition by Greiner et al. j2] the fleld has developed rapidly. 

Apart from the perfectly uniform lattice potentials generated by a single optical 
standing- wave, other lattice topologies have been discussed. Through the use of 
speckle patterns, optical lattices with random disorder have been realized j3]. Already 
the superposition of two standing-wave lattices of different wavelengths leads to a 
superlattice with a spatial modulation of the lattice well depths. For all of these non- 
uniform lattice topologies, rich phase diagrams have been predicted jSl El IZ|, which 
include new phases like an Anderson localized or a Bose-glass phase. 

The technique of sympathetic cooling of fermionic atoms through thermal contact 
with a Bose gas has paved the way to experiments with Bose-Fermi mixtures in optical 
lattices. The phase diagram of these systems of mixed quantum statistics shows a 
multitude of new phenomena |H1 El ^3 First experimental results ^2] using an 
admixture of a fermionic species to a Bose- Einstein condensate in an optical lattice reveal 
flrst effects of the inter-species interactions on the matter-wave interference pattern. In 
addition to the lattice topology, Feshbach resonances can be employed to tune the 
parameters of the Bose-Fermi-Hubbard Hamiltonian and thus to fully explore the zero- 
temperature phase diagram. 

The aim of this paper is twofold. First, we investigate the effect of non-uniform 
lattice potentials — a simple two-colour superlattice in our case — on the phase diagram 
of binary Bose-Fermi mixtures. Second, we introduce a physically motivated basis 
truncation scheme to extend the calculations to experimentally relevant system sizes. 
The basic formulation of the model and benchmark results for purely bosonic systems 
are summarized in sections |21 and El respectively. We use a direct numerical approach 
to solve the large-scale eigenproblem of the Bose-Fermi-Hubbard Hamiltonian, which 
constrains us to relatively small systems. In order to treat larger systems by direct 
diagonalization, we introduce an adaptive basis truncation scheme, which reduces the 
many-body basis to the physically relevant states. The truncation scheme is discussed 
in section |3] and a flrst application to bosons in superlattices is presented in section |S] 
Finally, the phase diagram of Bose-Fermi-mixtures in superlattices is discussed in detail 
in sections ini and [71 
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2. Bose-Fermi Hubbard Hamiltonian with a superlattice potential 

We consider a system of Na bosonic and Nc fermionic atoms in a one-dimensional 
optical lattice with / lattice sites at zero temperature. A detailed investigation of 
these systems in regular lattice potentials can be found in Refs. jHl In extension 
of these studies, we are going to discuss the impact of spatial lattice modulations on 
the phase diagram of boson-fermion mixtures. In the simplest case we assume a two- 
colour superlattice potential, generated by a superposition of two standing-wave lattices 
of different wavelengths. Provided the lattice is sufficiently deep, i.e., the gap between 
the ground and first excited band is sufficiently large, we can restrict ourselves to the 
lowest energy band jlSj. A straightforward generalization of the single-band Bose- 
Hubbard Hamiltonian leads to the single-band Bose-Fermi-Hubbard Hamiltonian. In 
second quantisation it reads [3 |H1 lEl HI] : 

H = - al,d, + clA + h.a.) + 1 Ka T^nf"^ {n^^^ - l) 

1=1 ^ 1=1 

+ V^^jzhf^nf + ^jZe, {hf+h^[^). (1) 
1=1 1=1 

The operators a] (a^) create (annihilate) a boson in the lowest Wannier state localized at 
lattice site The operators cj (ci) are the corresponding fermion creation (annihilation) 
operators, all fulfilling their particular commutation or anticommutation relations. The 
occupation number operators for the bosonic and fermionic component are given by 
n["'' = ajai and n^'^'' = cjci, respectively. 

The first term of the Hamiltonian (P) describes the tunnelling of atoms between 
adjacent lattice sites with a tunnelling matrix element J. We assume periodic boundary 
conditions, thus tunnelling between sites / and 1 is included. The next two terms 
describe the boson-boson and the boson-fermion on-site interaction with strengths Vaa 
and Vac, respectively. Due to the single-band description and the Pauli-principle there 
is no on-site fermion-fermion interaction. The last term describes an on-site potential 
energy in which A is the maximum amplitude and the e; specify the superlattice 
topology. The parameters of the Hubbard Hamiltonian encapsulate the properties 
of the many-body system, such as the depth and shape of the lattice potential and 
the scattering lengths of the boson-boson and the boson-fermion interaction. They 
are defined via matrix elements in the localized Wannier states. The presence of a 
superlattice potential is refiected mainly in the site-dependent single-particle energies 
ei. The site-dependence of J, Vaa and Vac can be neglected in good approximation. 

The many-body states within the Hubbard model can be conveniently represented 
in an occupation number basis. For a mixed boson-fermion system, the full occupation 
number basis is constructed via a direct product of the bosonic and the fermionic 
occupation number bases. A general many-body state can then be written as a 
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superposition of all possible combinations of the number states: 

I ^ ) = E E I {n'fK ... }. ) ® I {n'CK ... , 4^}, ) . (2) 

a=l /3=1 

The sets of occupation numbers {ni^\ ■■■ for the bosonic component comprise all 

possible distributions of the A^^^ atoms over the / lattice sites with = ^a- ^oi 

the fermionic occupation numbers {ni^\ ... , n^f^} the additional constraint n\'^^ G {0, 1} 
reflects the Pauli principle. The basis dimension grows rapidly with increasing lattice 
size and particle numbers, the dimensions of the bosonic and the fermionic number bases 
are given by Da = {Na + I- l)!/(Arj(J - 1)!) and = I\/{N^{I - iVj!), respectively, 
with the total dimension being the product of the two sub-space dimensions. 

Because of the factorial growth of the basis dimension, any exact solution of the 
many-body problem using the full basis is restricted to moderate system sizes. To solve 
the eigenvalue problem of the Bose-Fermi-Hubbard Hamiltonian, we are using refined 
Lanczos-type algorithms fHl, which enable us to treat basis dimensions of the order of 
10^. 

One of the major advantages of exact diagonalization techniques is their 
applicability throughout the whole phase diagram, from the coherent superfiuid regime 
to the strongly correlated Mott-insulator [21 E] . The price for this universality is the 
restriction to moderate system sizes due to the huge basis dimensions. Our aim is to 
increase the computationally feasible system sizes by a physically motivated, adaptive 
truncation scheme, which selects the important basis states for a certain region of the 
phase diagram out of the complete basis. This allows us to approach the experimentally 
relevant system sizes, to check for finite size effects, and to use the truncated basis in 
time-dependent simulations 

3. Single component Bose gas — a reminder 

We start with a brief reminder of some previous results [Zj on Bose gases in superlattices, 
which later on will serve as a benchmark for the accuracy of truncated calculations. 
Consider a single-component Bose gas with A^^^ = 10 particles in a two-colour superlattice 
with / = 10 lattice sites. The topology of the superlattice is defined by the set of on- 
site energies ei depicted in figure ^ The Hamiltonian follows from by omitting all 
fermionic contributions. Expressing all energies in units of the tunnelling energy J, 
the interaction strength Vaa and the superlattice amplitude A are the remaining two 
parameters. 

A first impression of the rich structure of the phase-diagram in the Vaa — A plane 
is provided by the maximum coefficient in the expansion of the groundstate: 

CL, = max{C^} (3) 

depicted in figure Efa). Note that the index j3 has been dropped due to the absence of 
the fermionic species. In the case of small values of C^^^ there are no preferred number 
states. For small superlattice amplitudes A and boson-boson interaction strengths Vaa, 
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Figure 2. (a) Maximum coefficient C^iax in the expansion of the number states (0) 
in the Vaa ~ A plane, (b) Maximum number fluctuations cimax- Both for Na = 10 
bosonic particles and / = 10 lattice sites. 



the system is in the superfluid phase, where virtually all number states contribute. 
If C^Qj. is close to 1, the groundstate is dominated by a single number state. In 
the Mott-insulator phase at large interaction strengths Vaa this is the number state 
with one particle per lattice site. Others are suppressed because of their interaction 
energy contribution. The homogenous Mott-insulating phase is present as long as the 
superlattice amplitude is smaller than the interaction energy A < Vaa as can be seen in 
figure Efa). As soon as they are approximately equal, A ^ Vaa, double occupations of 
the deepest lattice wells become energetically possible. In this transition region many 
number states contribute, giving rise to the reduced C^^^ at A ^ Vaa- This is the onset 
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of the quasi Bose-glass phase, where the interplay of lattice irregularities and interaction 
causes redistributions of the particles within the system. These redistributions of the 
particles also manifest in increasing number fluctuations 

(^liax = max {{4j\n^ - {4j\ni\tpf} (4) 

shown in figure Efb). Throughout the phase diagram one can clearly identify the 
complementary behaviour of C^^^ and the number fluctuations amax- The lobes result 
from the finite differences of the on-site energies in the superlattice potential (compare 
figure P). In the limit of an infinite lattice with random disorder these stable regions will 
vanish and the genuine Bose-glass phase will emerge. Another characteristic regime is 
the localized phase in which all particles are gathered at the deepest lattice well of each 
super-cell. It is obvious that this phase only exists at very small interaction strengths. 
The effect of localization due to disorder in the lattice is called Anderson localization 

m- 

4. Basis truncation scheme 

One possible way to extend these direct solutions of the eigenvalue problem to larger 
systems consists in a controlled and physically motivated truncation of the underlying 
basis. The challenge is to decide which number states are important and which can 
be discarded from the outset. As an example consider a system deep within the Mott- 
insulating phase. Number states with several particles at one particular lattice site are 
associated with very large interaction energies and their contribution to the ground state 
is negligible. Those less important number states can be omitted from the beginning, 
leading to a significant reduction of the basis dimension. 

A simple but very efficient a priori measure for the importance of individual number 
states is the expectation value of the Hamiltonian. The following inequality is used to 
select the relevant basis states: 

( a{ni, ...,ni} \ H \ {ni, nj}^ ) < Etmnc- (5) 

Hence, number states associated with interaction or on-site energy contributions above 
a specific truncation energy Etmnc are discarded. By varying the parameters Vaa, Vac, 
and A one can optimize the basis for a particular region of the phase diagram. By 
varying the truncation energy Etmnc one can adjust the size of the remaining basis and 
at the same time assess the truncation errors. The truncation scheme works particularly 
well in regions where the ground state is composed of a few dominant basis states, i.e. 
in regions of large C^^^. On the other hand, in regions where almost all number states 
contribute significantly, e.g. in the superfiuid phase, the truncation scheme will be less 
reliable. 

The generation of the truncated basis starts from a reference number state that 
minimizes the energy expectation value. With respect to this reference state, multiple 
particle- hole excitations are generated until the new states violate the inequality (jSJ. If 
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system 


sections (A/ J) 


dimension 


reference number state 


strong truncation 

PS 0.3% of complete basis 


1-28 


91 


1 {1,1,1,1,1,1,1,1,1,1} ) 


29-50 


288 


1 {1,1,2,1,0,1,1,2,1,0} ) 


51-150 


199 


1 {0,2,2,1,0,0,2,2,1,0} ) 


moderate truncation 
PS 4% of complete basis 


1-29 


3743 


1 {1,1,1,1,1,1,1,1,1,1} ) 


30-50 


1866 


1 {1,1,2,1,0,1,1,2,1,0} ) 


51-150 


2737 


1 {0,2,2,1,0,0,2,2,1,0} ) 


complete basis 


1-150 


92378 





Table 1. Basis dimension and reference states for different truncation levels and 
different ranges of A/ J. 



the truncation energy Etmnc is large enough, this process generates the complete set of 
number states. 

As a first benchmark, we study the effect of the basis truncation for the single- 
component Bose gas discussed earlier. For the sake of clarity and simplicity the 
interaction energy is fixed at Vaa/ J = 30. The calculations performed with the 
complete basis are compared to calculations with two different truncation levels. The 
phase diagram is split into three parts, each of them is calculated with an optimized 
basis. Tabled lists the different ranges in A/ J together with the corresponding basis 
dimensions and reference states. 

Since this approach fulfils the variational principle, we allow for overlaps at the 
borders of different sections and select the basis that yields the lowest groundstate 
energy. The strong truncation includes 3-particle-3-hole excitations with respect to an 
appropriate reference state chosen for the particular region in the phase diagram. The 
remaining dimension is about 0.3% of the complete basis. For the basis with moderate 
truncation, higher-order excitations are also included. Nevertheless, the dimensions of 
the moderate truncated basis is still reduced to 4% of the complete basis. 

The dependence of the maximum coefficient C^^x ^^e amplitude of the 
superlattice A at fixed interaction strength is depicted in figure a). If the amplitude 
A is smaller than the interaction strength Vaa, the system is in the homogenous Mott- 
insulator phase and the number state with one particle per lattice site dominates. 
Further increase of the superlattice amplitude leads to the onset of the quasi Bose- 
glass phase, where a pattern of regions with small values of C^^^ and large number 
fluctuations and regions with small number fluctuations and large values of C^^^ emerges 
as depicted in panels (a) and (b) of figure El The differences between the calculations 
with the full and the truncated basis are barely visible. Only for the strongly truncated 
basis with 0.3% of the full basis dimension, noticeable deviations occur. The basis using 
a moderate truncation to 4% of the full dimension yields results which are practically 
identical to the full calculation. 

A well-defined observable is the fraction of particles that undergo Bose-Einstein 
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Figure 3. Na = 10 bosonic particles on / = 10 lattice sites at fixed interaction 
strength Vaa/J ~ 30. (a) maximum coefficient C^^xj (b) maximum number fluctuation 
<^max, (c) condensate fraction fcond, a-U plotted over the amplitude of the superlattice 
A/ J. Calculations within the complete basis (solid line), truncated to 4% of the 
complete basis (dashed-dottcd, blue online), and truncated to 0.3% of the complete 
basis (dashed, red online). The results of the full basis and the truncated basis with 4 
% of the full dimension coincide almost everywhere. 



condensation. The largest eigenvalue Ai of the one-body density-matrix defined by: 

Pw = ( I 4«« I ) (6) 
is used to calculate the condensate fraction [71 QE] ■ 

fcond — ITT- ('^) 

The eigenvectors of the one-body density-matrix are called natural orbitals — in a 
regular lattice these are Bloch functions. 

Figure Efc) illustrates the dependence of the condensate fraction fcond on the 
superlattice amplitude A for fixed Vaa/ J = 30. Leaving the homogenous Mott-insulator 
phase, the condensate fraction increases rapidly, reaching a local maximum at the 
transition point to the quasi Bose-glass phase at A/ J = 30. After a slight decrease 
within the first lobe of the Bose-glass phase it increases further. It seems that in the 
limit of large superlattice amplitudes A, two separate "localized condensates" [12] form 
at the deepest lattice wells. The term condensate has to be used with caution in this 
context - we will come back to this point in section [HI 

Again, the deviations between truncated and full calculations are generally small. 
At the local maximum at A/ J = 30 the exact calculation yields fcond = 0.278, for the 
moderate truncation fcond = 0.276 and for the strong truncation fcond = 0.258. It is an 
intrinsic feature of the truncation to underestimate the condensate fraction because of 
the number states that are missed during the calculation of the one-body density-matrix 
elements. The kink in figure [Hfc) at A/ J = 50 (dashed line) is due to the change of the 
basis at this point (see table Q). In principle one could adopt another optimized basis 
in this region in order to improve the description. 

An observable that is directly accessible to the experiment is the matter-wave 
interference pattern after release from the lattice and free expansion. After switching 




Figure 4. Matter-wave interference pattern in the vicinity of the Mott-insulator to the 
quasi Bose- glass transition point. Calculations within the complete basis (solid line), 
truncated to 4% of the complete basis (dashed-dotted, blue online), and truncated 
to 0.3% of the complete basis (dashed, red onhne). The dashed-dotted hue is hardly 
visible because of the good reproduction of the complete calculation. 



off the confining potential, the intensity of the matter-wave at a point x can be written 

as m 

J(x) = A^{x)A{x) IV)- (8) 

Since we are not interested in the spatial envelope of the interference pattern, the 
amplitude operator A{x) depends only on the phase 0;(x) acquired on the way from 
site I to the observation point x, 

i(x) = 4^Ee^'^'(^)a,. (9) 

V J 1=1 

Considering microscopic distances between the lattice sites and a macroscopic distance 
from the lattice to the observation point, we assume a constant phase difference between 
adjacent sites 6 = 0/+i(x) — (f)i{x). This far-field limit leads to the following expression 
for the matter-wave intensity: 

m = \ E ( ^0 I I V^o ) = y E e^('-'')^ p^) , (10) 

i,i'=i i,i'=i 

in which we used equation ©. This expression is closely related to the quasi- momentum 
distribution. If 6 is equal to one of the possible quasi-momenta = ^irk/I, where 
k = 0,1,...,/ — 1, we obtain the respective quasi-momentum occupation numbers. 
Exemplarily figure |3] shows the matter- wave interference pattern for three values of the 
superlattice amplitude A in the vicinity of the transition point from the homogenous 
Mott-insulator to the quasi Bose-glass phase. From the sequence of interference patterns 
one can see that the height of the peak at quasi-momentum g = reaches a maximum at 
A = Vaa with an absolute value smaller than the condensate occupation number. Note 
that in a regular lattice the condensate wave function — the lowest natural orbital — is 
identical to the quasi-momentum zero Bloch function. Hence, in that case the occupation 
of the quasi-momentum zero mode is equal to the lowest eigenvalue of the one-body 
density matrix. This relation does not hold for irregular lattice potentials. 
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Figure 5. Left: Na = / = 15 system with 3 super-cells, right: Na = / = 20 
system with 4 super-cells, (a) and (b): Matter- wave interference pattern T{S) at the 
homogenous Mott-insulator to quasi Bose-glass transition point Vaa/J = A/J = 30. 
(c) and (d): condensate fraction fcond over the superlattice amplitude A. 



As already mentioned, the truncation scheme is expected to deviate most from the 
complete calculations if a large number of basis states is necessary to represent the 
ground state. But even in the transition region, where this occurs, less than 1% of the 
complete set of number states is sufficient to reproduce the properties of the system 
qualitatively. The basis including 4% of the complete set of number states included 
shows virtually no deviation from the exact calculation. It seems that the estimate 
(jni) is a suitable a priori measure to select the physically relevant number states. It 
is obvious that the more irregular the optical lattice is, the less energetically similar 
number states exist and the better the truncation scheme works. 

5. Multiple super-cells and finite size effects 

As a first application of the truncated basis in a regime where calculations with the full 
basis are not feasible anymore, we study the impact of finite size effects for a single- 
component Bose gas in a superlattice. To this end, we extend the two-colour lattice to 
three and four super-cells, corresponding to a problem of 15 bosons on 15 lattice sites 
{Na = I = 15) and 20 bosons on 20 lattice sites {Na = I = 20), respectively. We employ 
truncations to 1% and 0.001% of the respective full basis dimensions. 

The resulting matter wave interference patterns for Vaa/ J = A/ J = 30, i.e., at the 
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transition from homogeneous Mott-insulator to the quasi Bose-glass phase, are depicted 
in the upper part of figure El They can be compared to the result of the full calculation 
for a system of two super-cells shown in figure lU The matter- wave interference patterns 
seem to be barely affected by the truncation and are of the same general shape for all 
system sizes. Only the increasing number of possible quasimomenta qk = 27t/I with 
/c = 1, 2, ... J — 1 is reflected in the substructure of the interference peaks. 

Figure El also demonstrates the dependence of the condensate fraction on A/ J at 
fixed Vaa/ J = 30 for the truncated three and four super-cell system. The comparison 
with the complete calculation for the two super-cell system displayed in figure |21 reveals 
that the condensate fraction is generally reduced. The reason for this is two-fold: First, 
the basis truncation leads to an underestimation of the condensate fraction as discussed 
in section IH Second, finite size effects are relevant for the condensate fraction. The 
simplest indication for this results from the existence of a lower bound for fcond which is 
proportional to 1/J. The largest eigenvalue Ai of the one-body density matrix ((7|), which 
defines the condensate fraction, has a lower bound which results from the normalization 
of the one-body density matrix: 



Hence, the smallest possible value for the condensate fraction is given by =1/1. For 
the genuine Mott-insulating phase in a regular lattice of increasing size the condensate 
fraction would go to zero like 1/J. In agreement with this intrinsic size dependence, the 
condensate fraction within the Mott-insulating phase shown in figure El tends to smaller 
values. 

From these investigations we can conclude, that for the major observables finite 
size effects are under control and that calculations performed in small systems are well 
suited for qualitative and even semi- quantitative predictions for larger systems. 

6. Binary Bose-Fermi Mixtures 

Recent developments in the sympathetic cooling of boson-fermion mixtures in optical 
lattices pave the way to experimental studies of dilute, degenerate gases with mixed 
quantum statistics ^2] • Various aspects of boson-fermion mixtures in uniform lattices 
have been discussed before |2m jU] . In this section, we study the effect of an additional 
superlattice potential on the phase diagram. 

The following calculations were performed for a system with / = 10 lattice sites 
and Na = Nc = 5 particles of each species using the complete basis of Da ■ D^ = 504504 
number states. The boson-boson interaction energy is fixed to Vaa/ J = 20 in order to 
simplify the discussion — other (non-zero) values would not change the basic structure 
of the phase diagram but only cause a rescaling. 

Figure El shows the maximum coefficient C^ax a function of the boson-fermion 
interaction strength Vac/ J and superlattice amplitude A/ J. Already the maximum 




i=l 
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Figure 6. Left panel: Cf-^^^^ of a binary Bose-Fermi mixture with Na = Nc ~ 5 
particles and / = 10 lattice sites. Right panel: The mean occupation numbers (white 
for bosons and gray for fermions) for the selected points in the phase diagram. The 
boson-boson interaction energy is fixed at Vaa/J = 20. 



coefficient of the number state representation of the ground state reveals the rich 
structure of the phase diagram. We can identify several straight lines or valleys of 
minimum C^^^ that subdivide the phase diagram into different regions. There is 
a horizontal line at A = Vaa and a vertical line at Vac = Vaa, whose locations are 
determined by the boson-boson interaction strength and which are independent of the 
precise actual lattice topology. In contrast, the additional diagonal valleys are governed 
by the details of the lattice topology, i.e., the set of on-site energies e^. 

Typical occupation number distributions for bosons and fermions in the different 
areas of the phase diagram separated by those valleys are also presented in figurelHl They 
correspond to the values of the boson-fermion interaction strength and the superlattice 
amplitude marked by the labels in the phase-diagram. The characteristics of the different 
regions can be summarized as follows: 

(A) In this region, the atoms are uniformly distributed over the lattice because of 
the dominant repulsive interactions. All number states with few particles per 
site are favoured since the strong boson-fermion repulsion suppresses states with 
different species occupying the same site. Because of the commensurate filling 
and the strongly repulsive interactions between bosons as well as between bosons 
and fermions this phase is closely related to the Mott-insulating phase in a purely 
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bosonic system. 

(B) The mean occupation numbers of both species follow the shape of the superlattice 
potential. In the upper-left corner of this region two bosons are able to share the 
deepest superlattice wells with one fermion. Increasing boson-fermion repulsion Vac 
forces one boson to leave the deepest superlattice well and the mean occupation 
numbers exhibit a profile as depicted in inset (B). 

(C) Above A/ J = 20 = Vaa/J the on-site energy of the deepest superlattice wells 
exceeds the boson-boson interaction energy and they are capable to hold 2 bosonic 
particles. This double occupation occurs if the boson-fermion interaction prevails 
the boson-boson interaction and forces the fermion to leave the deepest well. 

(D) Complete localization of the bosons at the deepest superlattice sites induced by the 
strong repulsive boson-fermion interaction despite of the boson-boson repulsion. In 
a purely bosonic system with Vaa/ J = 20, a full localization would occur only for 
extremely strong superlattice potentials with amplitudes A/J > 210. The presence 
of the fermions and the strong boson-fermion repulsion, causing a separation of 
bosons and fermions, promote the localization. The borders of this region of 
separation induced localization depend not only on the interaction strengths and 
the superlattice amplitude, but also on the filling. The more particles per super-cell 
are available, the further the onset of localization is pushed to larger superlattice 
amplitudes. 

(E) Transition region between (B) and (C). This area of subtle rearrangements is 
governed by the finite differences of the on-site energies ei of the superlattice 
potential. 

(F) Transition region between (B) and (D). The fermions on sites / = 2 and 7 
successively force the bosons to the state of complete localization. If the bosons 
were the only species this would happen at much larger superlattice amplitudes. 

The phase diagram in terms of the condensate fraction fcond of the bosonic species 
is depicted in figure [7|^a). Only three areas with a significant condensate fraction 
are visible. In the lower left corner — at small boson-fermion interaction strengths 
and superlattice amplitudes — large condensate fractions occur despite of the strongly 
repulsive boson-boson interaction {Vaa/J = 20), since the incommensurate boson filling 
averts a Mott-insulator phase. The condensate is depleted as soon as either the 
increasing boson-fermion interaction leads to a Mott-like phase, or the superlattice 
tends to localize the particles. The large condensate fractions in an extended region 
in the upper right part and along the upper part of the Vac/ J = axis have to be 
treated with care and are of the same character as the localized condensates discussed 
in section |3] The almost complete localization of the bosonic particles at the deepest 
lattice sites / = 3 and 8 manifests in the corresponding diagonal matrix elements, i.e., 
the occupation numbers, of the one-body density matrix pgg"* and pgg^ In comparison, 
the other matrix elements are small except for the off-diagonal elements p^^^ and p^^^ . 
This results in two large eigenvalues of the one-body density matrix and, therefore, in 
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Figure 7. Binary Bose-Fermi mixture with Na ~ Nc ~ 5 particles on / = 10 lattice 
sites at fixed boson-boson interaction energy Vaa/J = 20. (a) condensate fraction 
fcond of the bosonic species, (b) fringe visibility ly of the bosonic species. Although the 
condensate fraction at the top right is large, the fringe visibility vanishes there. This 
is due to the different nature of the localized condensates (see text). 



our description, in a large condensate fraction. Note, however, that this 'condensate' 
is completely different from the condensate associated with the superfluid phase in a 
homogenous lattice, where all particles are delocalized and exhibit long-range phase 
coherence. This characterization of the condensate is closely related to the existence of 
off-diagonal long-range order ^H] in the one-body density matrix, which is absent in the 
case of the localized condensates. 

This interpretation is also supported by the behaviour of the fringe visibility ly 
shown in figure Hl^b). The visibility is directly derived from the matter- wave interference 
pattern pUj) and defined as: 

-'-max -'-mm /in\ 
-'-max \ -'^mm 

Clearly, the visibility is large at small boson-fermion interaction strengths and 
superlattice amplitudes (lower left corner) as expected for a conventional condensate. In 
contrast to this, there is no sizable fringe visibility in the upper-right region, indicating 
that the localized phase is incoherent and has to be distinguished from the genuine 
Bose-Einstein condensate. This also applies to the narrow stripe along the Vab/ J = 
axis, where the fringe visibility exhibits an intermediate value. Again the groundstate is 
dominated by number states that have many bosonic particles at the deepest superlattice 
wells. In contrast to the area at the top right, not all of the particles are localized and 



Bose-Fermi mixtures in ID optical superlattices 



15 



thus many more number states contribute. The more number states contribute to the 
groundstate, the more the phase coherence is restored leading to more pronounced 
interference fringes. 

In a recent experiment, Giinter et ah jT2] found that the successive admixture of 
fermions to a Bose-Einstein condensate in an optical lattice leads to a reduction of the 
interference fringes of the bosonic component. A detailed inspection of the visibility of 
the interference fringes u in figure [7|^b) at A/J = might confirm their observation. 
With increasing interaction strength Vac the visibility also decreases rapidly. Albeit the 
boson-fermion interaction strength is fixed and the ratio between the particle numbers, 
Na/Nc, is changed in the experiment, the effect on the bosons should be similar to the 
setup presented here. 

7. Basis truncation for Bose-Fermi mixtures 

The basis truncation scheme that was introduced in section |3J can also be applied to 
Bose-Fermi mixtures. Before going to larger systems, we repeat the calculations of the 
previous section for a I = 10, Na = Nc = 5 boson-fermion mixture using a truncated 
basis. The truncation of the basis again makes use of the importance criterion (jSj). For 
the present example the basis dimension is reduced to Da - Dc = 72762, thus about 15% 
of the number states of the complete basis are kept. We use this rather large set because 
the phase diagram is not partitioned into several areas with different optimized bases 
— for simplicity, only one basis optimized for A = Vaa = Vac is employed. Of course, 
one can also work with an optimized basis for each point or for particular regions of the 
phase diagram. 

For brevity, we only consider the condensate fraction and the fringe visibility, as 
shown in figure |Hl for the truncated calculations. The general structure of the phase 
diagram of figure [7| is well-reproduced with the truncated basis. The absolute value of 
the condensate fraction is underestimated in the truncated calculation as was discussed 
earlier. In close connection to this, the narrow stripe at the top left of figure IHl^b) shows a 
significantly smaller visibility u than in the complete calculation. This is because number 
states with many particles at the same lattice sites are discarded by the truncation 
scheme due to the optimization to a region, where also the boson-fermion interaction 
plays a significant role. Although the coefficients Ca/s of the missing number states 
would have been small, their contribution to the off-diagonal elements of the one-body 
density matrix is sizable. Their omission causes a reduction of coherence and thus the 
reduced visibility of the interference fringes. 

Despite the quantitative deviations in some special regions in the phase-diagram, 
its global structure can be fully described by the truncated basis. Naturally, the use of 
adapted bases in the different characteristic regions in the phase diagram would improve 
the quantitative agreement. 
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Figure 8. Binary Bose-Fermi mixture with Na = Nc = h and / = 10 for fixed 
Vaa/ J = 20 calculated with a truncated basis of roughly 15% of the full dimension, 
(a) condensate fraction fcond of the bosonic species, (b) fringe visibility v of the 
bosonic species. The comparison with figure [3 illustrates the quality of the truncated 
calculation. 

8. Bose-Fermi mixtures beyond half- filling 

Finally, we investigate a system with Na = Nc = 7 particles and J = 10 lattice sites. 
The complete basis for this system consists of Da ■ Dc = 1372800 number states. For 
the presented calculations a truncation to 22% of the complete basis is used. 

The maximum coefficient C^^^, depicted in figure IHl as a function of Vac/ J and 
A/ J for fixed Vaa/ J = 20, reveals a phase diagram with a structure different from 
the half-filling case shown in figure IHl The most significant difference is the absence 
of the Mott-insulator-like phase in the lower part of the phase diagram. Even in 
the limit of very strong repulsive interactions, a homogeneous Mott-insulator phase 
is not possible, since there are always lattice sites occupied with more than one 
particle due to the incommensurate filling. For this reason in region (A), where the 
boson-fermion interaction is negligible compared to the boson-boson repulsion, pairs 
of bosons and fermions occupy the lattice with a distribution following the topology 
of the superlattice. An increase of the superlattice amplitude (region (B)) leads to a 
bosonic double occupancy of the deepest superlattice wells, whereas the fermions are 
forced to spread over several sites with mean occupation numbers less than or equal 
to one. The horizontal line, which appeared in the half-filling case for superlattice 
amplitudes A equal to the boson-boson interaction strength Vaa, vanishes because of 
the incommensurability. For strong boson-fermion repulsion, any slight irregularity 
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Figure 9. Boson-fcrmion mixture with / = 10, Na = Nc = 7. Left panel: maximum 
coefficient C^^j^. in the expansion of number states, right panel: mean occupation 
numbers nf" ('^■1'^^) of the bosonic (fermionic) species. The boson-boson interaction 
energy is fixed at Vaa/ J — 20. 



of the lattice — i.e. a small nonzero superlattice amplitude — leads to multiple boson 
occupancies at the deeper lattice sites. A comparison of the mean occupation numbers 
at (C) and (D) or at (E) and (F) shows that if the boson-fermion interaction exceeds 
the boson-boson interaction, the mean occupation numbers are stable against increasing 
superlattice amplitude. However, at the same time the increasing value of C^^^ indicates 
that the number fluctuations decrease successively. 

9. Conclusions 

We have studied the phase diagram of binary boson-fermion mixtures in optical 
superlattices. For the symmetric half-filled system we identified a homogenous Mott- 
insulator-like phase at small superlattice amplitudes and strong repulsive interactions 
within the bosonic species and between the bosonic and the fermionic component. A 
phase of complete localization of the bosonic species can be realized already for moderate 
superlattice amplitudes if the boson-fermion repulsion clearly exceeds the boson-boson 
repulsion. These "localized condensates" are of a different nature than the genuine 
Bose-Einstein condensates in optical lattices as indicated by the absence of matter-wave 
interference fringes. 

The phase-diagram of a boson-fermion mixture beyond half-filling is less structured, 
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since the homogeneous Mott-insulator phase at small superlattice amplitudes is absent. 
In the region of strong boson-fermion interaction any slight irregularity within the 
lattice immediately leads to double occupancies of the deeper site. For increasing 
superlattice amplitude, the mean occupation numbers are almost unaltered, only the 
number fluctuations decrease. 

Part of the calculations were performed using a controlled basis truncation scheme 
tailored to select the physically important basis states out of the complete set. This 
versatile tool allows to increase the particle numbers and lattice sizes accessible to 
diagonalization techniques. By explicit comparison with full calculations, we have shown 
that truncated bases containing only about one percent of the number states, already 
provide quantitatively reliable results. The simplicity of the basis truncation facilitates 
time-dependent calculations for systems of moderate size, as will be discussed in a 
separate publication [TH] . 
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